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Abstract. The aim of this introductory article is two-fold. First, we aim to offer 
a general introduction to the theme of Bose-Einstein condensates, and briefly discuss 
the evolution of a number of relevant research directions during the last two decades. 
Second, we introduce and present the articles that appear in this Special Volume of 
Romanian Reports in Physics celebrating the conclusion of the second decade since the 
experimental creation of Bose-Einstein condensation in ultracold gases of alkali-metal 
atoms. 


1. INTRODUCTION 

1.1. ATOMIC BOSE-EINSTEIN CONDENSATES 

The Bose-Einstein condensate (BEC) is a macroscopic quantum state of matter, 
which was predicted theoretically by Bose and Einstein 90 years ago [1]. Atomic 
BECs were created experimentally in ultracold vapors of s, Rb [2], 23 Na [3] and 
'Li [4] 70 years later. The aim of this Special Issue is to celebrate the twentieth 
anniversary of this remarkable achievement, which was also marked by the Nobel 
Prize in Physics for 2001, awarded to E. A. Cornell, W. Ketterle, and C. E. Weiman 

[5]. 

In the BEC state, all atoms in the bosonic gas fall (“condense”) into a single 
quantum-mechanical ground state. The transition to the BEC occurs if the atomic 
density, n, and the de Broglie wavelength. A, corresponding to the characteristic 
velocity of the thermal motion of the atoms, satisfy the following condition [6, 7]: 

nA 3 > 2.612, 


( 1 ) 
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which implies that A is comparable to or larger than the mean distance between 
atoms, thus making the gas a macroscopic (degenerate) quantum state. In the above- 
mentioned atomic gases, this theoretical condition is met at temperatures T which 
are a small fraction of milli-Kelvin (mK), hence the atomic BECs created in the lab¬ 
oratories are, as a matter of fact, the coldest objects existing in the universe (the early 
BEC-experiments achieving the condensed state around 100 nK). Their creation be¬ 
came possible after the development of appropriate experimental techniques needed 
to reach the necessary ultra-low temperatures (see, e.g., Ref. [8]). The required ex¬ 
treme cooling is achieved in two stages. First, the method of laser cooling (which 
was also rewarded with the Nobel Prize in Physics for 1997 [9]) is applied to the 
gas loaded into a magneto-optical trap, which makes it possible to create a moder¬ 
ately cool state, at temperature ~ 100 fiK. Next, this state undergoes forced evapora¬ 
tive cooling, losing ~ 90% of atoms, and the remaining atomic cloud spontaneously 
forms the BEC. In the experiments, the number of atoms in the BEC typically ranges 
between 1,500 and 1,000,000 (although both smaller and larger numbers are, in prin¬ 
ciple, possible), and the size of the domain in which the gas is trapped is ~ 100 //m. 
A characteristic time scale relevant to the experiments is measured in milliseconds, 
while the lifetime of the condensate can be easily raised to several seconds. 

The applicability of the laser-cooling method to particular atomic species de¬ 
pends on the peculiarities of their electron configuration. As a result, this technique 
has made it possible to achieve the Bose-Einstein condensation in vapors of alkaline, 
alkaline-earth, and lanthanoid metals: 'Li, 23 Na, 39 K, 41 K, 85 Rb, 87 Rb, 133 Cs, 52 Cr, 
40 Ca, 84 Sr, 86 Sr, 88 Sr, 174 Yb, 164 Dy, and 168 Er. Perhaps especially interesting, among 
the more recent developments, is the creation of BEC in the gases of chromium [10] 
and dysprosium [11], where atoms carry large magnetic moments, which makes it 
possible to predict and observe many effects produced by long-range dipole-dipole 
interactions [12]. The challenging aim of creating BEC in the gas of spin-polarized 
hydrogen atoms has been finally achieved too, with a specially devised technique 
which made it possible to cool the gas to 50 pK [13]. 

The BECs created in the laboratory constitute a prototypical manifestly quan¬ 
tum-macroscopic state of matter available in the experiments. In other settings, where 
low temperatures are crucially important too, macroscopic quantum effects, such as 
superconductivity in metals and superfluidity in liquid helium, are well known too, 
but they correspond to “implicit” quantum states. For instance, a superconducting 
metallic sample as a whole is not a macroscopic quantum object. The same pertains 
to the recently created out-of-equilibrium BEC of quasi-particles in condensed mat¬ 
ter, namely, exciton-polaritons [14] and magnons [15], which have drawn a great deal 
of attention in the past decade (see Refs. [16-18]). Also, a considerable attention was 
drawn to the topic of localization of exciton-polaritons in semiconductor microcavi¬ 
ties [19-22]; for an excellent recent review focused on several physical phenomena 
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exhibited by exciton-polariton condensates see Ref. [23]. The condensation of effec¬ 
tively massive photons trapped in a microcavity was reported too [24], the peculiarity 
of these settings being the nonconservation of the total number of the quasi-particles 
or photons. 

Surveys of the broad subject of Bose-Einstein condensation and numerous re¬ 
lated areas are provided by many review articles and books [6, 7], [12], [25]-[41]. It 
is important to mention that this list of surveys on the topic of BEC is, of course, far 
from being exhaustive. This is a clear indication of the impact of this research theme 
to almost all branches of contemporary physics. 

The goal of the present article is to offer a broad picture of some of the past and 
currently active research areas in the realm of BEC (admittedly biased towards the 
particular research interests of the authors) and to overview the scientific literature, 
akin to a Resource Letter of American Journal of Physics. 


1.2. MEAN-FIELD DESCRIPTION AND NONLINEAR DYNAMICS OF BEC 

From a theoretical standpoint, and for many experimentally relevant conditions, 
the static and dynamical properties of a BEC can be described by means of an effec¬ 
tive mean-field equation known as the Gross-Pitaevskii (GP) equation [6, 7]. This 
is a variant of the famous nonlinear Schrodinger equation (NLSE), incorporating an 
external potential used to confine the condensate; NLSE is known to be a universal 
model describing the evolution of complex field envelopes in nonlinear dispersive 
media [42-44]. In the case of BECs, the nonlinearity in the GP (NLSE) model is 
introduced by the interatomic interactions, accounted for through an effective mean 
field. Thus, an inherent feature of the BEC dynamics is its nonlinearity, which is in¬ 
duced by collisions between atoms, in spite of the fact that the density of the quantum 
bosonic gases is very low. 

The studies of the matter waves in the presence of the nonlinearity drive a 
vast research area known as “nonlinear atom optics” (see, e.g., Refs. [45]). Impor¬ 
tantly, many collective excitations, including self-trapped localized states supported 
by the condensate’s intrinsic nonlinearity (e.g., solitons), are less straightforward to 
create (and difficult to describe by adequate models) in dense media featuring macro¬ 
scopic quantum phenomena, such as liquid helium. An exception are phase solitons 
(fluxons), i.e., quanta of magnetic flux trapped in long Josephson junctions formed 
by superconductors [46], whose experimental and theoretical studies are relatively 
straightforward and have been developed in detail [47]. Nevertheless, atomic BECs 
constitute an ideal setting for studies of such macroscopic nonlinear excitations, as is 
explained in more detail below. 

In the atomic BEC with intrinsic self-attraction (e.g., 7 Li or 85 Rb BECs), the 
creation of effectively one-dimensional (ID) matter-wave bright solitons (both iso- 
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lated ones and multi-soliton sets) in cigar-shaped configurations, which are tightly 
confined by external potentials in the transverse plane, was successfully reported in 
condensates of 'Li [48-50] and 85 Rb [51, 52] (see also the reviews in Ref. [53]). 
Collisions between moving quasi-ID solitons also admit accurate experimental im¬ 
plementation [54] and theoretical analysis [54, 55]. 

More typical is the repulsive sign of the inter-atomic interactions (as in the 
cases of s ' Rb or 23 Na BECs), which lends the BEC the effective self-repulsive non¬ 
linearity. This kind of the intrinsic interaction readily creates dark solitons, which 
were predicted theoretically [56] and created experimentally [57, 58] in BECs loaded 
into a cigar-shaped trap. In fact, dark solitons were first created [57] prior to the re¬ 
alization of the above-mentioned bright solitons in the self-attractive condensates, 
placed into the same type of the trapping potential (a review of the topic of dark 
solitons in BEC was given in article [36]). Similar to the case of bright solitons, 
not only single-soliton [57] states, but also multiple dark solitons were created [58], 
while their interactions and collisions were also studied both in theory [58, 59] and 
experiments [58, 60]. 

Later, stable dark-bright soliton complexes in binary BEC were predicted in 
theory [61] and observed in experiments [62] as well. In more recent works, multiple 
dark-bright solitons [63], as well as dark-dark solitons [64] were also experimentally 
created. In addition, in the same setting of the self-repulsive nonlinearity, not only 
dark solitons, but also bright solitons are possible too: in particular, if - instead of 
the usual parabolic trap - a periodic (optical lattice) potential [31] is used to confine 
the condensate, then gap solitons can be formed. Experimentally, matter-wave gap 
solitons built of ~ 250 atoms in a 8 'Rb condensate, were reported in Ref. [65]. 

On the other hand, in the case of multidimensional BEC geometry, there has 
been an intense theoretical and experimental activity on vortices and vortex structures 
in BECs with the self-repulsive intrinsic nonlinearity (see Refs. [66] for reviews on 
this topic). This is due to the fact that vortices are intimately related to the superfluid 
properties of BECs, and play an important role in transport, dissipative dynamics 
and quantum turbulence (see, e.g., Refs. [67]). Historically, the first observation of 
vortices in BECs was achieved by phase imprinting [68] (a technique that was also 
used to create dark solitons [57]). Nevertheless, there exist other techniques that 
have been used in experiments to nucleate vortices in BECs: these include stirring 
the condensate above a certain critical angular speed [69] (this method was used 
to create vortex lattices [70]), nonlinear interference between different condensate 
fragments [71] (this technique was also employed for the creation of dark solitons 
[58]), by forcing superfluid flow around a repulsive Gaussian obstacle within the 
BEC [72], and through the Kibble-Zurek mechanism [73] (the latter was originally 
proposed for the formation of large-scale structures in the universe by means of a 
quench through a phase transition [74]). 
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It is also important to mention that there still exist other notable results con¬ 
cerning nonlinear phenomena in BECs. These include the realization of few-vortex 
clusters and complex structures such as vortex dipoles, vortex tripoles, parallel vor¬ 
tex rings, etc., and the study of their dynamics [69, 72, 75-88], the observation of 
quantum shock waves [89], the realization and study of “hybrid” soli ton-vortex struc¬ 
tures [90], the observation of Josephson oscillations [91] and spontaneous symmetry 
breaking transitions [92] in BECs loaded into a symmetric double-well trapping po¬ 
tential, and so on. 

Many of the nonlinear phenomena mentioned above can be successfully studied 
in the framework of mean-field theory. Nevertheless, there exist other phenomena 
(where the underlying intrinsic BEC nonlinearity is also important) that can not be 
described by means of the mean-field theory. Examples of such situations, along with 
cases where quantum and/or thermal fluctuations come into play, will be discussed 
below. 


1.3. WHY ARE BOSE-EINSTEIN CONDENSATES ALWAYS “IN VOGUE”? 

A great advantage offered by the BEC in low-density atomic gases is that these 
media can be easily and very efficiently controlled by means of external magnetic 
and optical fields. This circumstance enables various experiments, and provides a 
framework for very accurate theoretical models. As a result, the ultracold gases 
can be used for emulation of many phenomena which originate, e.g., in condensed 
matter physics, but take a very complex form in the original settings, due to the strong 
interactions in them, while in atomic BECs similar effects can be simulated in a much 
simpler (“clean”) form. 

An example which has recently drawn a great deal of interest is the emulation of 
the spin-orbit coupling (SOC), i.e., the interaction between the motion of an electron 
or hole and their spin, or, in other words, the linear mixing between two components 
of the electron's or hole’s spinor wave function. It is a fundamentally important 
phenomenon in semiconductors, known in two distinct forms, as the Dresselhaus 
and Rashba SOC [93]. In binary atomic BECs, created as a mixture of two different 
hyperfine states of atoms of 8 ‘Rb, the SOC has been experimentally implemented 
as a similar linear coupling between the two atomic components, with the help of 
appropriate laser beams illuminating the condensate, and a dc magnetic field applied 
to it [94]. 

There are many other prominent examples of such an emulation provided by the 
atomic-gas BEC. One such example concerns the direct observation of the transition 
between the superfluid and Mott-insulator states formed in the BEC loaded into an 
“egg-carton” periodic potential for individual atoms, which is, in turn, imposed by an 
optical-lattice structure [95]. Other examples include resonantly-enhanced tunneling 
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in periodic potentials [96], artificial gauge fields [41, 97], topological insulators [41], 
the fractional quantum Hall effect [97], and so on. Recently, the general topic of 
using the BEC in atomic gases as the “quantum simulator”, and the related theme of 
using ultra-cold fermionic gases for the same general purpose, have been reviewed 
in a series of articles [98] and in a recent book [39], 

As concerns applications of BEC, the most straightforward one is, arguably, 
the use of matter waves for interferometry. With the help of the intrinsic nonlin¬ 
earity, they may feature superb accuracy, in comparison with traditional optical in¬ 
terferometric schemes [99]. Especially promising are interferometric schemes of 
the Mach-Zehnder type, based on splitting and subsequent recombination of bright 
matter-wave soli tons [53, 100, 101]. More generally, nonlinear waves that arise in 
BECs play a principal role in different applications, in addition to their intrinsic in¬ 
terest. For instance, bright solitons have been argued to provide the potential for 
100-fold improved sensitivity for interferometers and their lifetime of a few seconds 
enables precise force sensing applications [102] (see also Ref. [103]). For repul¬ 
sively interacting BECs too, the inter-atomic interactions have also been suggested 
to increase the sensitivity to phase shifts, precisely due to the emergence of dark 
solitons which enable (e.g. through their oscillation in the confining trap) better de¬ 
tection schemes [104]. Furthermore, nonlinear atom interferometers can overcome 
the limitations of the current state-of-the-art standard based on the so-called Ram¬ 
sey spectroscopy [105, 106], due to their ability to surpass the classical precision 
limit. Finally, vortices present their own potential for applications. An example 
is the so-called “analogue gravity” [107], whereby they may play a role similar to 
spinning black holes. This allows to observe, in experimentally controllable envi¬ 
ronments, associated phenomena such as the celebrated Hawking radiation or even 
simpler ones such as the super-radiant amplification of sonic waves scattered from 
black holes [108]. 

BECs have also been experimentally demonstrated to be usable for perform¬ 
ing remarkable tasks such as the implementation of the Gaussian sum algorithm for 
factoring numbers [109], by exploiting higher order quantum momentum states, im¬ 
proving in this way the algorithm's accuracy, once again, beyond its classical imple¬ 
mentation. This is in line with the development of the Shor algorithm as an efficient 
quantum mechanical way to factorize large numbers, a task thought to be classically 
intractable [110, 111]. 

Another potentially promising application is the use of BEC as a resource 
for the implementation of quantum computing. In this context, one possibility is 
to use trapped droplets of the condensate as qubits [112, 113]. In optical lattices 
also, atomic analogs of semiconductor electronic circuits (the so-called “atomtron- 
ics”) have been proposed, in order to realize quantum devices such as diodes and 
transistors [114]. On the other hand, collision between quantum matter-wave soli- 
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tons can be used to induce entanglement between them [115], which is a prerequisite 
necessary for the design of soliton-based quantum-information-processing schemes. 

Also promising is the development of atom-wave lasers, which should be able 
to emit high-intensity coherent matter-wave beams, in continuous-wave (CW) or 
pulsed (soliton-like) regimes. Such beams may be very useful, in particular, for pre¬ 
cision measurements. The first experimental realization of a CW matter-wave laser 
was reported in Ref. [116], which was followed by the development of a design with 
a separated BEC reservoir and the beam-emitting cavity [117]. These experimental 
works used condensates consisting of 87 Rb atoms. A review of experimental results 
on the topic of matter-wave lasers was recently published in Ref. [118]. Models of 
matter-wave lasers operating in a pulsed regime were developed theoretically [119]. 


2. MODELS, SETTINGS AND BASIC RESULTS 


2.1. THE GROSS-PITAEVSKII EQUATION 


As mentioned above, the fundamental model which provides for an accurate 
description of the BEC in dilute degenerate gases of bosonic atoms is based on 
lowest-order mean-field theory. According to this approach, the gas is described 
by means of the Gross-Pitaevskii equation (GPE) for the single-particle wave func¬ 
tion, T ( x,y,z,t ), where ( x,y,z) and t are the coordinates and time [as mentioned 
above, “degenerate” means that the de Broglie wavelength of atoms moving with 
the thermal velocity in the dilute gas is large enough in comparison with the mean 
inter-atomic distance - see Eq. (1)]. The three-dimensional (3D) form of the GPE is 
written, in physical units, as: 



) ^ + U(x 1 y,z-t)'i + ^-a s (x,y,z-t)\^\ 2 ^, 
J m 

( 2 ) 


where m is the atomic mass, U (. x,y,z;t ) is the external potential acting on indi¬ 
vidual atoms and thus confining the condensate as a whole (U may depend on time 
too, which is often called management of the potential [120]), and a s (x,y,z;t ) is 
the scattering length which determines collisions between the atoms: a s > 0 and 
a s < 0 correspond to repulsive and attractive interactions, respectively. The spatial 
and temporal dependence of a s , which is important for many predictions of non¬ 
trivial dynamical states in the BEC (see below), may be induced by means of the 
Feslibach-resonance (FR) management technique [120]. FR implies the formation 
of a quasi-bound state of two atoms in the course of the collision between them in 
the presence of an external magnetic field [121], or under appropriate laser illumina¬ 
tion [122]. The FR can be induced too by combined magneto-optical settings [123]. 
Making use of the very accurate tunability of the FR [124], spatially non-uniform 
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and/or temporally variable external fields, controlling the FR, can be employed to 
induce spatially and/or time-dependent nonlinearity coefficients, accounted for by 
a s (x,y,z;t). Further, the wave function is subject to the normalization condition, 
which is determined by the total number, N, of atoms in the condensate: 



| 'F (x, y, z; t) | 2 dxdydz = N. 


( 3 ) 


Note that, alternatively, the wave function may be defined with the unitary norm, 
f11 1^ {x,y,z-,t)\ 2 dxdydz = 1, replacing a s in Eq. (2) by Na s . 

The GPE and its variants constitute a relatively simple mathematical frame¬ 
work, which admits precise simulations and the use of effective analytical approaches. 
The latter include the Thomas-Fermi approximation (TFA), which neglects the terms 
of the second derivative (the kinetic energy of the quantum particles) in Eq. (2) 
[7], and a more accurate and versatile variational approximation, which has found 
a great number of applications to BECs [125, 126]. The TFA is relevant for a s >0 
(the self-repulsive nonlinearity), for solutions which, in the simplest case, do not 
include a nontrivial phase structure; in such a case, the solution with chemical po¬ 
tential y > 0 is found by using the ansatz T = e _ ' / ' / < I > 'f j7,\ (x. y. z), where the density 
I^tfa^,//;")! 2 is approximated as: 


|$TFA0,y,2)| 2 


m 

4irh 2 a s ( x,y,z ) 


fJ-~U(x,y,z), for U(x,y,z) < y, . 
0, for U(x,y,z) > y. 


The TFA may be easily generalized for vortex states, which arc sought for as solu¬ 
tions to Eq. (2) in the cylindrical coordinates (p = x 1 + y' 2 . 8, z), assuming that 
U = U (p, z) and a s = a s ( p, z) are subject to the cylindrical symmetry: 

T = e - lflt+lSe ® (p, z ), (5) 


with real integer vorticity S, and function <t> satisfying the respective stationary equa¬ 
tion: 


y$ = - 


2m 


d 2 


d 2 1 d__S^_ _ 

dp 2 ^ p dp p 2 dz 2 


<f> + U(p,z)$ + ^-a s (p,z) |<f>| 2 <f«. (6) 


m 


The TFA neglects all the terms with p- and ^-derivatives in Eq. (6), which yields the 
following generalization of solution (4) [127, 128]: 


|^tfa| 2 = 


m 


47t h 2 a s ( p,z ) 


( /i ” 

U(.p,z) + <^-p~ 2 

l 

0, for U (p, 


2 m 


2 m 


( 7 ) 

Techniques for numerical treatment of GPEs have been developed in great de¬ 
tail too. They include methods for the solution of boundary-value problems, aimed at 
finding stationary states trapped in external potentials, or self-trapped due to the non¬ 
linearity (solitons), as well as direct simulations of the GPE in real or imaginary time 
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(the latter approach helps to generate stationary solutions for ground states [129], due 
to the relaxational character of the dynamics). In particular, the semi-implicit split- 
step Crank-Nicolson algorithm [130] has become a method of choice for solving 
GPEs in many settings. The unconditional stability of this algorithm makes it espe¬ 
cially useful in studies of real-time dynamics of BECs, although it can be equally 
well used to produce ground states of relevant experimental setups, including fast- 
rotating BECs with many vortices. The readily available Fortran and C codes [131], 
which implement the Crank-Nicolson approach in ID, 2D, and 3D geometries of 
BECs with different symmetries, are well tested and highly optimized. Furthermore, 
the C codes are parallelized using the OpenMP approach, which speeds up execu¬ 
tion of numerical simulations of BECs significantly, up to one order of magnitude on 
modern computers with multi-core CPUs. Apart from the imaginary-time propaga¬ 
tion implemented in the framework of GPEs, ground states of a BEC (and of other 
quantum systems described by linear and nonlinear Schrodinger equations) can also 
be calculated numerically using the higher-order effective-action approach [132], as 
demonstrated in Refs. [133]. Other popular methods used to solve the GPE by means 
of pseudospectral and finite-difference methods are detailed in Refs. [134-137]. A 
survey of analytical and numerical methods used in the studies of BEC models was 
given in review article [32], 

For the BEC composed of chromium [10] or dysprosium [11] atoms, with mag¬ 
netic moments, /a, polarized in certain direction by an external uniform dc magnetic 
field, the GPE includes a nonlinear term accounting for the long-range dipole-dipole 
interactions: 





1 — 3 cos 2 0 


dr'. (8) 


r — r 


In Eq. (8), V 2 stands for the 3D Laplacian, the term ~ a s , which represents the usual 
contact interaction (see Eq. (2)), is dropped on grounds that it is small in comparison 
to the dipole-dipole interaction, //q is the magnetic permeability of vacuum, and 6 is 
the angle between vector r — r' and the polarization direction of the atomic magnetic 
moments. 

Beyond the mean-field approximation, quantum fluctuations and interaction 
of the condensate with the thermal component of the gas are described within the 
framework of the Hartree-Fock-Bogoliubov equations, which are essentially more 
cumbersome than the relatively simple GPE [138]. However, it is only under some¬ 
what special conditions (i.e., small atom numbers below N ~ 1000 or large enough 
temperatures, of the order of many tens or hundreds nK) that fluctuations play a 
crucial role for coherent matter-wave patterns, including solitons in the case of ex¬ 
periments. Nevertheless, such settings are becoming of increasing interest in both 
theoretical and experimental studies [37, 40]. 
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2.2. TWO-COMPONENT SYSTEMS 

Shortly after the experimental realization of the single-component BEC, ad¬ 
vances in trapping techniques opened the possibility to simultaneously confine atomic 
clouds in different hyperfine spin states. The first such experiment, the so-called 
pseudospinor condensate, was achieved in mixtures of two magnetically trapped hy¬ 
perfine states of 87 Rb [139]. Subsequently, experiments in optically trapped 23 Na 
[140] were able to produce multicomponent condensates for different Zeeman sub- 
levels of the same hyperfine level, the so-called spinor condensates. In addition to 
these two classes of experiments, mixtures of two different species of condensates 
have also been created by sympathetic cooling (i.e. condensing one species and al¬ 
lowing the other one to condense by taking advantage of the coupling with the first 
species) - see examples of such BEC mixtures below. 

Regarding the modelling of multi-component BECs, it is natural to proceed 
from the single-component GPE to the corresponding system of coupled GPEs. In 
particular, for the simplest case of two-component mixtures, Eq. (2) is replaced by 
the following system of equations for mean-field wave functions, T' 1 and Tg, of the 
two components: 

V 2 ^ i + C/ i ( r ;t)^ i + ^-(4 1) |^ i l 2 + 4 12) |^2| 2 )^t, (9) 

ot Ami mi ' / 

V 2 ^2 + U 2 (r;t) tf 2 + (a^ |^ 2 | 2 + 4 12) |^i| 2 ) ^ 2 , (10) 

Ot 27712 m 2 V ' 

( 12 ) 

where a s is the scattering length for collisions between atoms belonging to the two 
different species, for which the trapping potentials, U\ and U 2 , induced by the same 
external field, may be, generally speaking, different. Masses m\ and m 2 are different 
for BEC mixtures composed by different atom species - also-called heteronuclear 
systems - such as 85 Rb— 8 'Rb [141] and s 'Rb— 133 Cs [142] binary BEC mixtures. 
On the other hand, m\ = m 2 for mixtures of different hyperfine states of the same 
atomic species; such mixtures were experimentally realized for the first time with 
8 'Rb atoms [143]. In the case of repulsive intra- and inter-species interactions, when 

all the scattering lengths in Eqs. (9) and (10) are positive, that is, of > 0, of > 0, 

( 12 ) 

and a s >0, the condition for the immiscibility of the two components - and the 
onset of the separation between them - is given by [144]: 

4 1} 4 2) < (°4 12) ) ’ (U) 

which implies that the repulsion between atoms belonging to the different compo¬ 
nents is stronger than the repulsion between atoms in each component. Condi¬ 
tion (11) pertains to the free infinite space, while the pressure of the trapping potential 
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makes the binary BEC more miscible, shifting the critical point, = 

( 12 ) 

to larger values of positive a s [145]. Immiscible two-component condensates, 



loaded into a trap, form domain walls separating the two components [146]. Such 
domain walls were observed in experiments [147], 

Numerical analyses of realistic trapped states of binary BECs have revealed, 
in the case of cigar-shaped geometries, two distinct immiscible stationary configura¬ 
tions: a segregated one, in which the two components face one another, being sep¬ 
arated by a domain wall, and a symbiotic state, in which one component effectively 
traps the other [148]. Both configurations do not directly obey the aforementioned 
simple miscibility criteria, although they can be transformed into a miscible configu¬ 
ration when the condensate is subject to a resonant drive [148]. Finally, we mention 
that, while the symbiotic state is not a soliton, numerous types of solitons are known 
to exist in spinor condensates (see, e.g., [149, 150]). 

In the case when the two wave-function components represent different hy- 
perfine states of the same atomic species, an external resonant radiofrequency field 
(with frequencies in the GHz range) may add linear mixing, with strength k (alias 
Rabi coupling), to the system, which is accounted for by extra terms and k x \> \, 
added to equations (9) and (10), respectively [151]. The interplay of the Rabi cou¬ 
pling with the repulsive interactions causes a shift of the miscibility-immiscibility 
transition (11) [152] (see relevant experimental results in Ref. [153]). 


2.3. TRAPPING POTENTIALS 


Coming back to the single GPE (2), it is relevant to stress that two most com¬ 
mon types of the confining potential are the harmonic-oscillator (HO) 



( 12 ) 


and optical-lattice (OL) 


U OL = U x sin 2 (k x x) + U y sin 2 (k y y) + U z sin 2 (k z z) 


(13) 


ones, where fl 2 2 are trapping frequencies of the (generally, anisotropic) HO po¬ 
tential, and U x , y ,z represent the depths of the periodic OL potential. The OL is built 
as the classical interference pattern by pairs of counterpropagating mutually coher¬ 
ent laser beams illuminating the condensate, with respective wavelengths \ x , y , z = 

2 t /k x ,y,z- 

The OL is made attractive or repelling for individual atoms by red- or blue¬ 
detuning of the illuminating light with respect to the frequency of the dipole transition 
in the atoms. Typically, the OL wavelength A ~ 1 pm is used in experiments. Usu¬ 
ally, the depth U of the OL potential is measured in natural units of the recoil energy, 
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Er = ( hk ) 2 / (2m). The use of the OL potentials for the creation of matter-wave pat¬ 
terns in BEC was proposed in [154] and relevant applications, such as macroscopic 
quantum interference, immediately ensued [155] - see also reviews [156] and [157]. 
A well-known example of the effect induced by the OL is the transition from the 
bosonic superfluid to the Mott insulator [95]. 

Generally, the technique based on the use of OLs is similar to that which was 
proposed [158] and implemented in the form of photonic lattices in photorefractive 
media, producing a number of spectacular results, including ID and 2D optical soli- 
tons and vortices in 2D [159] - see reviews [160]. 

2.4. THE DISCRETE SYSTEM 

In both the BEC and photonic settings, a very deep (compared to the chemical 
potential) OL potential effectively splits the mean-field matter-wave function, or the 
optical electromagnetic field, into a set of nodes (each one representing one well) 
weakly interacting between them via tunneling coupling. In this case, using the ex¬ 
pansion of the continuum field over a set of Wannier modes localized around local 
wells, GPE (2) can be reduced to a discrete nonlinear Schrodinger equation (NLSE) 
[161, 162]. In a properly scaled form, its 3D version is 

. d '& j i. i 1 

i—= ~ ^ {^j+l,k,l + ^j-l,k,l + ^j,k+l,l + ^k,k-l,l 

+ ^j,k,l+l + ^j,k,l-l — ^j,k,l)— \^fj,k,l\ 2 ^j,k,h ( 14 ) 

where j,k,l are discrete coordinates on the lattice, and x l' h kj are amplitudes of 
trapped matter-wave fragments, the 2D and ID versions being produced by obvi¬ 
ous reductions of Eq. (14). The present form of the discrete NLSE implies that 
the on-site nonlinearity is self-attractive. However, unlike the continuous model, 
in the discrete one the self-repulsive nonlinearity may be transformed into its self- 
attractive counterpart by means of the well-known staggering transformation [162], 
^j,k,l{t) = (—1 y +k+l f r l2l/ T* k ft), where the asterisk stands for the complex con¬ 
jugate. 

The discrete NLSE gives rise to many species of soli tons [162], especially 
interesting ones being discrete localized vortices, which were predicted theoreti¬ 
cally [163] and created experimentally (as nearly discrete objects) in photonics using 
waveguide arrays built in a photorefractive material [159]. Thus, while ID and 2D 
versions of Eq. (14) apply to the photonic settings [160], the 3D discrete system is 
meaningful solely in the BEC context. In the latter case, it generates complex sta¬ 
ble localized modes, such as, e.g., discrete Skyrmions [164], diamonds, octupoles, 
oblique vortices, and vortex cubes [165], among many others. 

The discrete NLSE suggests a direct transition to the fully quantum system 
in the form of the Bose-Hubbard (BH) model, replacing the mean-field (classical) 
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lattice wave functions in Eq. (14) by quantum operators, bj (for simplicity, discussed 
here in the ID setting). The corresponding Hamiltonian is 


H = Y1 ~ Jh ) ( 6 7+i + 6 i-i) + - Un i ( n 7 - !) 


(15) 


where rij = //■ bj is the operator of the on-site number of particles, J is the inter- 
site-hopping constant, and U is the constant of the on-site interaction (U > 0 and 
U < 0 correspond, as before, to the self-attraction and self-repulsion, respectively). 
A famous result produced by Hamiltonian (15) is the phase diagram separating the 
quantum superfluid and Mott insulator (see Ref. [166]). In terms of applications, 
the BH model is a natural tool for the theoretical analysis of operations of BEC- 
based qubits. Reviews of the topic of BH in connection to its realization in BEC 
can be found in articles [34] and [113]. The well-elaborated numerical technique for 
the analysis of the BH model and its modifications is based on the density-matrix- 
renormalization-group method [167]. 


2.5. REDUCTION TO LOWER-DIMENSIONAL SETTINGS 

The above-mentioned nearly-ID cigar-shaped traps are represented by the po¬ 
tentials which include tight confinement in the transverse plane, i.e., large f l 2 = H 2 = 
Q 2 _ in Eq. (12), and an arbitrary weak potential, U(x,t), acting in the axial direction, 
x. In this case, the 3D wave function can be approximated by the factorized Ansatz 
(see, e.g., Ref. [168]), 

1 ( y 2 T z‘ 2 \ 

®(x,y,z,t) = -j=^— exp - —p—J (16) 

where the transverse part represents the ground state of the two-dimensional HO in 
the transverse plane, with the respective oscillator length 

a± = s/hf (_ l ) (17) 

(typical values relevant to the experiments are a± ~ 3 //m), while the axial wave 
function, tp{x,t), subject to the normalization condidon f 1 x \ip(x)\ 2 dx = N, see 
Eq. (3), satisfies the ID equation obtained by averaging in the transverse plane: 

M h 2 d 2 i!> . . 4h 2 , .. l2 

m m + - x ' tm 

which has the form of the one-dimensional NLSE, with an external potential, U. 
Essentially the same equation occurs in many other physical settings, such as the 
light propagation in planar waveguides, in which case t is actually the propagation 
distance, while —U{x) represents the confining prodle of the local refractive index 
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[43, 44, 169]. Accordingly, Eq. (18) with a s < 0 and a s > 0 gives rise, respectively, 
to the commonly known bright- and dark-soliton solutions. Interestingly, the NLSE 
in ID is integrable in the case of U = 0 and a s = const [170]. 

An interesting ramification of this setting is the toroidal quasi-ID trap, which 
is described by Eq. (18) with periodic boundary conditions in x. Such toroidal traps, 
realized by means of several different techniques, are available in the experiment 
[171]. 

A quasi-2D pancake-shaped (oblate) configuration, with strong confinement 
acting in the transverse ID direction, corresponds to large Q'j, = fE in Eq. (12), 
combined with a general relatively weak potential, U (x,y), acting in the pancake’s 
plane. This configuration is approximated by the respective factorized ansatz, 

1 f i z 2 \ 

y{x,y,z,t) = ^ 1/4 ^ / _ exp (19) 

which leads to the effective 2D equation, 



2m 


f d 2 d 2 \ 
\ dx 2 dy 2 ) 




2V2nti 2 

ma _1 


a s (x;t) IV’I 2 '*/’- 


( 20 ) 


The reduction of the 3D GPE (2) to its ID version (18) on the basis of factorized 
Ansatz (16) with the fixed transverse localization radius, a±, is relevant in the limit of 
low density. For higher density, the reduction is also based on ansatz (18), in which 
ax is allowed to be a variable parameter. Then, the reduction to ID is performed by 
means of the variational approximation, which leads to a system of ID equations for 
i/j (x,t) and a± (x) [172], that can be reduced to a single effective equation for the 
ID wave function with a nonpolynomial nonlinearity. The resulting “nonpolynomial 
NLSE” (abbreviated as NPSE [172]), and the respective local expression for a±, are 
given (in a scaled form) by 


and 



1 d 2 ^ 

2 dx 2 


f U(x,t)ipA 


l + ( 3 / 2),#| 2 , 

xA+s'lVi 2 


( 21 ) 


a\ = 1 + glV’l 4 , g = 2a s /a± < 0. (22) 

However, the relevant reduction from 3D to ID may be done in multiple ways (e.g., 
working at the level of underlying Lagrangian/Hamiltonian structure and of the cor¬ 
responding action or at that of the equations of motion). On a related direction, using 
the standard adiabatic approximation and accurate results for the local chemical po¬ 
tential, one obtains an alternative equation with a nonpolynomial nonlinearity [173]: 



ld 2 ^ 


f U (x,t)ip + 



2 dx 2 


(23) 
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where g is the same as in Eq. (22). 


2.6. COLLAPSE OF ATTRACTIVE CONDENSATES 


In the free space (U = 0), with a constant negative scattering length, corre¬ 
sponding to the self-attractive nonlinearity (a s < 0), a rescaled form of Eq. (20) 
amounts to the 2D version of the NLSE: 



1 / d 2 

2 \<9x 2 



1>- IV’I 2 ^- 


(24) 


A well-known fact is that Eq. (24) gives rise to a family of isotropic, so-called 
Townes’ solitons [174], 

ip = exp (— igt) (p(r), r = s/x 2 + y 2 , (25) 


with arbitrary chemical potential g < 0, and real function (p obeying the equation: 

g4> + (1/2) (4>" + r -1 <//) +4> 3 = 0. (26) 

The family of the Townes’ solitons is degenerate, in the sense that their norm takes a 
single value which does not depend on g: 

/»oo 

N T = 2n cp 2 (r)rdr «5.85. (27) 

Jo 

Note that an analytical variational approximation for this norm predicts Nt = 2ir, 
with relative error ~ 7% [175]. 

On the other hand, the three-dimensional GPE (2) in the free space, with the 
uniform self-attractive nonlinearity, (7 = 0 and a s < 0, gives rise to a family of 
isotropic soliton solutions in the form given by Eq. (25). Unlike their 2D counterparts 
in the form of the above-mentioned Townes’ solitons, the norm of the 3D solitons de¬ 
pends on g, N = const - (— g)~ 1 ^ 2 , cf. Eq. (27). The celebrated Vakhitov-Kolokolov 
(VK) necessary stability criterion [176, 177], dN/dg < 0, does not hold for this 
N(j.i) dependence, hence the entire family of the 3D free-space solitons is unstable, 
which is completely corroborated by the full analysis of the stability [177]. For the 
2D Townes solitons, Eq. (27) formally predicts neutral VK stability, dN/dg = 0, but 
in reality the Townes solitons are unstable too. However, their instability is nonlinear, 
i.e., it is not accounted form by any unstable eigenvalue in the spectrum of eigen- 
modes computed around the stationary soliton, using the respective Bogoliubov-de 
Gennes equations (BdGEs). The eigenvalue associated with the instability in this 
special case is at 0, corresponding to a special invariance arising in this critical case, 
namely the so-called conformal invariance [43], which allows a rescaling of the soli¬ 
tary wave. In fact, the instability of the multidimensional solitons is explained by 
the fact that the NLSE with the self-attractive cubic nonlinearity gives rise to the 
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dynamical collapse, i.e., the self-similar formation of a true singularity after a finite 
evolution time. In the 2D space, the collapse is critical, which implies, inter alia, that 
the norm of collapsing solutions must exceed a threshold (minimum) value, which 
is precisely the Townes-soliton norm (27), while the 3D collapse is supercritical, as 
its threshold norm is zero [177]. In the experiments with the self-attractive BEC, the 
onset of the collapse was readily observed (the first time in 85 Rb [178]), as sponta¬ 
neous explosion of the condensate (the so-called “Bose nova”). It is interesting to 
mention that the small part of the condensate surviving the explosion, can form a 
stable soli ton in 85 Rb [51]. 

2.7. MODELS FOR NON-BEC ULTRACOLD GASES 

The above discussion pertains to quantum bosonic gases. Ultracold fermionic 
gases have also been created in the experiment [179], which was followed by the ob¬ 
servation of their condensation into the bosonic gas of Bardeen-Cooper-Schrieffer 
(BCS) pairs [180]. The theoretical description of the Fermi gases is more com¬ 
plex, because, in the general case, the Pauli principle prevents the application of 
the mean-field approximation to fermions, making it necessary to treat such gases 
directly as systems with many degrees of freedom of individual particles (see, e.g., 
the works [181] and the review [182]). 

An approach to sufficiently dense Fermi gases is possible in terms of a hy¬ 
drodynamic description, which, in a sense, is a variety of the mean-field theory. This 
approach starts from the famous works by Yang and Lee who derived the energy den¬ 
sity for a weakly coupled BCS superfluid [183]. In the spirit of the hydrodynamic 
approach, an effective equation for an order parameter of the Fermi gas, T (x,y,x,t) 
was derived, which seems as the NLSE with the self-repulsive term of power 7/3 
[184]: 

f)tti fr 2 h 2 

* 6*7 = -~- V 2 * + U(x,y,z)* + — |*| 4 / 3 4/, (28) 

ot 2m e g 2m 

where m e g is the effective mass, which may be different from the particle’s mass. 
This equation is valid for a slowly varying order-parameter field, under the condition 
that the local Fermi energy is much larger than all other local energy scales, such 
as potential U ( x,y,z ). Equation (28) and its ID and 2D reductions can be used 
to predict various stationary and quasi-stationary density patterns in the Fermi gas 
[184], Related equations for the dynamics of Fermi gases are discussed in Ref. [185]. 

Similarly to the case of bosonic gases, nonlinear excitations of Fermi gases 
have attracted attention. In particular, dark solitons in a Fermi gas were predicted 
near the BEC-BCS transition, using the description in terms of the BdGEs [186], 
while their dynamical properties were studied in several works [187]. It is worth 
noting that dark solitons [188] and hybrid soliton-vortex structures [189] (the so- 
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called solitonic vortices) were recently observed in experiments with Fermi gases in 
the BEC-BCS crossover. 

Another well-known example of a dilute quantum medium different from BEC 
is the ID Tonks-Girardeau (TG) gas, formed by hard-core bosons, a solution for 
which may be mapped into that for a gas of non-interacting fermions [190]. In par¬ 
ticular, this mapping (also known as “Bose-Fermi mapping”) makes it possible to 
produce a solution for a dark soliton in the TG gas [191]. Importantly, the TG gas 
was realized experimentally using ultracold 87 Rb atoms loaded into a tightly confined 
quasi-ID hap [192], 

An equation for the order parameter of the TG gas was derived in work [193], 
in the form of the ID NLSE with the self-repulsive quintic term: 


in—— = — 


h 2 d 2 ^> 


+ U(x)^ + 


irti 2 


T'| 4 'F. 


(29) 


dt 2m dx 2 ' v ' ' 2m 

Similar to Eq. (28), Eq. (29) is valid only as a quasi-stationary one and it does 
not provide correct description of dynamical effects in the TG gas [194]. There 
exist examples of a relevant use of Eq. (29), including the prediction of solitons 
supported by the long-range dipole-dipole attraction between atoms forming the gas 
[195], and the study of dark soliton oscillations [196]. Interestingly, the result for the 
soliton oscillation frequency, which was analytically found in Ref. [196] to be equal 
to the axial hap frequency, was in agreement with numerical predictions obtained in 
Ref. [197] via the Bose-Fermi mapping. 


3. SOME SPECIAL TOPICS THAT HAVE ATTRACTED INTEREST 
3.1. TEMPORAL AND SPATIAL MANAGEMENT 

The management concept can be applied for the trapping potentials by making 
the HO or OL strengths time-dependent. A typical example of results produced by 
the time-periodic management of the HO potential (12), with O 2 = Oq + O 2 sin (ut) 
(in the simplest ID setting), is the prediction of parametric resonances of self-trapped 
matter-wave packets (solitons) in the latter setting [198]. On the other hand, for the 
OL the management was experimentally realized [199] in the form of a “rocking” 
OL, by introducing a small wavelength mismatch between the two laser beams build¬ 
ing the OL, which is made a periodic function of time: AA = AAosin(wi), i.e., the 
lattice as a whole performs periodic oscillatory motion. In particular, the rocking OL 
potential may effectively suppress the matter-wave tunneling across the lattice [199]. 

One particular application of the periodic modulation of the strength of the HO 
potential deals with the appearance of regular patterns in the density profile of the 
condensate through a modulational instability. A prototypical example is the emer¬ 
gence of Faraday patterns in cigar-shaped BECs [200] through the periodic modu- 
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lation of the strength of the radial component of the magnetic trap, similar experi¬ 
mental results being known for 4 He [201]. As another aspect of theoretical results, 
we mention numerous studies of Faraday waves in models of the condensates with 
short-range interactions [202-204], dipolar condensates [205], binary condensates 
with short-range interactions [206], Fermi-Bose mixtures [207], and superfluid Fermi 
gases [208]. As a related topic, let us mention that it has been shown theoretically 
that Faraday waves can be suppressed in condensates subject either to resonant para¬ 
metric modulations [209] or space- and time-modulated potentials [210, 211], and 
that pattern-forming modulational instabilities lead to chaotic density profiles [202] 
akin to those of turbulent BECs [67, 212]. Apart from the use of parametric excita¬ 
tions, the formation of density waves has been studied in expanding ultra-cold Bose 
gases (either fully [213] or only partly condensed [214, 215]), and the spontaneous 
formation of density waves has been reported for antiferromagnetic BECs [216]. 

As predicted theoretically (see Ref. [217] for a review), many possibilities for 
the creation of matter-wave patterns in BEC are offered by various patterns of spatial 
[218] and temporal [219] modulation of the local scattering length, a s , as implied 
by the general form of GPE (2). Experimentally, spatial “landscapes” of the scat¬ 
tering length can be induced, via the FR mechanism, by the corresponding spatially 
periodic distributions of the magnetic field (magnetic lattices), created with the help 
of periodic structures built of ferromagnetic materials [220]. Another possibility is 
the local modulation of the scattering length, which may be imposed, via the optical 
FR, by time-average patterns “painted” by rapidly moving laser beams [221], Spatial 
modulation of interatomic interactions has also been demonstrated at the submicron 
level via pulsed optical standing waves in an ytterbium BEC [222]. Once again in 
this context, it is relevant to point out that some of these possibilities, such as the 
temporal modulation of the nonlinearity have been also realized in parallel, in other 
contexts such as, e.g., nonlinear optics [223]. 

Employing such magnetically or optically induced Feshbach resonances, via 
the above-mentioned FR-management technique, indeed constitutes one of the most 
promising methods for manipulating BECs. Such a possibility to control the effective 
nonlinearity of the condensate, has given rise to many theoretical and experimental 
studies. Probably the most well-known example between these, is the formation of 
bright matter-wave solitons and soliton trains in attractive condensates [48]-[52], by 
switching the interatomic interactions from repulsive to attractive. 

This inspired many theoretical works studying the BEC dynamics under tem¬ 
poral and/or spatial modulation of the nonlinearity. In particular, the application of 
FR-management technique, with the low-frequency modulation of the strength of the 
magnetic field causing the nonlinearity to periodically switch between attraction and 
repulsion, can be used to stabilize 2D solitons in the free space [219] (in reality, a 
weak two-dimensional HO trapping potential is necessary in the experiment). How- 
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ever, this method does not work in 3D, nor for 2D vortex solitons. On the other hand, 
the same technique can be used for the generation of robust matter-wave breathers 
[224], 

On the other hand, the so-called “collisionally inhomogeneous condensates” 
(a term coined in Refs. [218]) controlled by the spatially modulated nonlinearity, 
have been predicted to support a variety of new phenomena. This new regime can 
be achieved by means of magnetically or optically controlled Feshbach resonances. 
The magnetic Feshbach resonances is a well-established experimental method, which 
has been used to study the formation of ultracold molecules [225], the BEC-BCS 
crossover [226], and the production of Efimov trimer states [227], but the inhomo¬ 
geneity length scale of the necessary magnetic field is usually larger than the size of 
the atomic sample, therefore this method is not very efficient in reaching the colli¬ 
sionally inhomogeneous regime. Flowever, the optical Feshbach resonance has been 
shown to allow fine spatial control of the scattering length, and recent experimental 
results demonstrate controllable modulations of the 5-wave scattering length on the 
scale of hundreds of nanometers [222]. 

The range of new nonlinear phenomena specific to the inhomogeneous-non- 
linearity regime includes adiabatic compression of matter waves [218, 228], Bloch 
oscillations of matter-wave solitons [218], emission of the solitons and design of 
atom-beam lasers [229], dynamical trapping of matter-wave solitons [230, 231] en¬ 
hancement of transmissivity of matter waves through barriers, [231, 232], creation of 
stable condensates exhibiting both attractive and repulsive interatomic interactions 
[230], competition between incommensurable linear and nonlinear lattices [233], the 
generation of dark solitons and vortex rings [234], control of Faraday waves [235], 
and many others. Importantly, the Feshbach resonance was used in recent exper¬ 
iments to induce real spatial inhomogeneities of the scattering length [221, 222], 
which paves the way for implementation of the above-mentioned phenomena in the 
experiment. 


3.2. MULTIDIMENSIONAL LOCALIZED STRUCTURES 

3.2.1. Attractive BEC 

The stabilization of multidimensional solitons is a problem of great interest not 
only to BEC, but also to nonlinear optics and related research areas [236-240]. It 
was predicted theoretically, but not as yet demonstrated experimentally, that the use 
of OL potentials is a universal means for the stabilization of such solitons [241, 242] 
(a similar stabilization mechanism was predicted for 2D optical solitons supported 
by the Kerr self-focusing nonlinearity in photonic-crystal fibers [243]). In particu¬ 
lar, the stabilization of 2D solitons is explained by the fact that the norm of such 
solitons, trapped in the OL potential, takes values below the threshold value (27) 
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corresponding to the Townes’ soliton, hence the solitons have no chance to start col¬ 
lapsing. The same property, N < At, explains the stabilization of solitons trapped 
in the HO potential (12) [244, 245]. For a comprehensive study of stability of two- 
dimensional elliptic vortices in self-attractive Bose-Einstein condensates, trapped by 
an anisotropic harmonic trapping potential, see Ref. [246]. 

On a different but related direction let us mention that OLs of the Bessel type 
can support soliton rotation, see Refs. [247-250] for the theoretical works and Ref. 
[251] for the experimental verification, a nonlinear phenomenon which has been in¬ 
vestigated in both BEC and nonlinear optical settings. 

Still more challenging objects are multidimensional vortex solitons (i.e., self- 
trapped modes with embedded vorticity). In addition to the possibility of the collapse, 
they are still more unstable to fragmentation by azimuthal perturbations [217, 236]. 
An accurate analysis of the stability of three-dimensional solitons with vorticity 
S = 1 in self-attractive Bose-Einstein condensates, trapped in an anisotropic three- 
dimensional HO potential was reported in Ref. [252]. The analysis predicts that 
vortex solitons can also be stabilized by OL potentials [241]. Of course, the lattice 
breaks the axial symmetry, but, nevertheless, the vorticity embedded into a localized 
state may be defined in this case too [241, 253]. In the simplest form, stable vortex 
solitons with topological charge (integer vorticity) S can be constructed as iV-peaked 
ring-shaped patterns with the vorticity represented by the phase circulation along the 
ring, with phase shift A (j> = 2ir/N between adjacent peaks (see the straightforward 
definition (5) of the vorticity for isotropic settings). OLs may stabilize solitons even 
if the lattice’s dimension is smaller by one than the dimension of the physical space, 
including ID [254] and 2D [254, 255] OLs in the 2D and 3D settings, respectively. 

A new approach to the creation of stable 2D solitons supported by the cubic 
self-attraction, which was considered impossible until recently, was put forward in 
theoretical work [256]. It is based on the system of two GPEs, linearly coupled by 
first-derivative terms representing the above-mentioned SOC of the Rashba type, and 
by nonlinear terms accounting for collisions between atoms belonging to the two 
different atomic states, which underlie the coupled system. In the scaled form, the 
system is given by 


and 


.dipi _ 

1 dt 


,d ^ 2 _ 

1 dt 


ivVi-(|^i | 2 + 7l^2| 2 )^i + A(^-z^) (30) 

^ 2 -(|^ 2| 2 + 7l^i|V2-A(^ + ^) , (31) 


where V 2 is the Laplacian acting on coordinates (x,y), the real-valued A is the 
strength of the SOC, and 7 is the relative strength of inter-component nonlinearity in 
comparison with the intra-component self-attraction. Due to the specific form of the 
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SOC terms, composite solitons are generated by the system (30)-(31) as bound states 
of a fundamental localized state in one mode and a vortex in the other mode ( semi- 
vortices , also referred to as filled vortices [257] or vortex-bright solitary waves [258]), 
or mixtures of fundamental and vortical components in both modes. As mentioned 
above, solitons trapped in the OL potential are stabilized by it because their norm 
drops below the threshold necessary for the onset of the collapse, while it was be¬ 
lieved that this is impossible in the 2D free space. The new feature of the composite 
solitons produced by the system (30)-(31) is that their norm also takes values below 
the threshold, without the help of any trapping potential. The full stability is provided 
by the comparison of values of the Hamiltonian of Eqs. (30)-(31), 



H 


T c.c. > dxdy , (32) 


where c.c., as well as the asterisk, stands for the complex conjugate, for the composite 
semi-vortex and mixed-mode solitons. The self-trapped modes of the former and 
latter types are stable and realize the system’s ground state (i.e., they minimize energy 
(32)) at 7 < 1 and 7 > 1, respectively. 


3.2.2. Repulsive BEC 


For the repulsive interaction between atoms, it has been predicted that stable 
matter-wave vortices are supported by condensates loaded into OLs [259] and that 
gap solitons and gap-soliton vortices also exist if the condensate is loaded into the 
OL of the same dimension [260]. On the other hand, the application of an OL, 
or of a magnetic lattice [ 220 ], to impose spatial modulation of the scattering length, 
a s (x,y,z) in Eq. (2), induces an effective nonlinear lattice, which can readily support 
ID solitons, while the stabilization of their 2D counterparts by nonlinear periodic 
potentials is a difficult problem [217]. 

New perspectives for the creation of stable complex 3D localized modes, such 
as vortex rings, vortex-antivortex hybrids, and Hopfions (twisted rings, featuring two 
independent topological numbers), which were unavailable in other physical media, 
were recently predicted by a model with the local strength of the self-repulsive cubic 
nonlinearity growing from the center to periphery at any rate faster than r 3 [128, 261, 
262], Realization of these settings in BEC is a challenge to the experiment. Below 
we will focus on some of the more standard settings involving vortices and related 
structures in higher-dimensional BECs without spatial or temporal modulation of the 
scattering length. 
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3.3. VORTICES AND VORTEX CLUSTERS 

Arguably, one of the most striking features of BECs is the possibility of sup¬ 
porting vortices, which have been observed in many experiments by means of a va¬ 
riety of methods. Vortices are characterized by their non-zero topological charge S, 
whereby the phase of the wavefunction has a phase jump of 2irS along a closed con¬ 
tour surrounding the core of the vortex. The width of single-charge vortices in BECs 
is of the order 0(£) - where £ is the healing length of the condensate (see, e.g., [7]) - 
while higher-charge vortices, with \S\ > 1, have cores wider than the healing length. 
Such higher-charge vortices are generally unstable in the homogeneous background 
case; nevertheless, they may be stabilized by employing external impurities [263] 
(the latter can be used to confine the so-called persistent current; see, e.g., the recent 
discussion of relevant experiments in [264] and references therein), or by using ex¬ 
ternal potentials [265]. Notice that, when unstable, higher-charge vortices typically 
split in multiple single-charge vortices, since the system has no other way to dispose 
of the topological charge [266]. 

The fact that single-charge vortices carry topological charge renders them ex¬ 
tremely robust objects: indeed, continuous deformation of the vortex profile cannot 
eliminate the 27r phase jump. An exception is a case where the background conden¬ 
sate density is close to zero, and that is why, in BEC stirring experiments, vortices 
are nucleated at the periphery of the harmonically trapped condensate [267]. 

Vortices are prone to motion caused by gradients in the density (and phase) of 
the background, induced by an external potential (as, e.g., in the case of a trapped 
BEC) or by the presence of other vortices. The motion of the vortex in such cases 
can be studied by means of the matched asymptotic expansion method [268]. The 
same method can also be used to study the effect of vortex precession induced by 
the external trap (see, e.g., the review [66]), also in the presence of collisional in¬ 
homogeneities and dissipative perturbations [269]. Note that in the simplest case of 
a single vortex in a 2D BEC confined in a harmonic trap, a Bogoliubov-de Gennes 
analysis reveals the connection of the vortex precession frequency with characteris¬ 
tic eigenfrequencies of the spectrum and - in particular - with the anomalous mode 
[270] (for the latter, the integral of the norm x energy product is negative [7]). In¬ 
deed, the negative-energy mode bifurcates in the linear limit from the dipole mode 
(which has a constant magnitude equal to the trap frequency); then, as the chemi¬ 
cal potential increases, the anomalous mode eigenfrequency decreases, and becomes 
equal to the precession frequency of the vortex in the Thomas-Fermi limit (see, e.g., 
Refs. [269, 271]). 

On the other hand, the motion induced on a vortex by another vortex is tan¬ 
tamount to the one observed in fluid vortices (see, e.g. Ref. [272]): this way, vor¬ 
tices with same charge travel parallel to each other at constant speed, while vortices 
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of opposite charges rotate about each other at constant angular speed. Following 
Helmholtz’ and Kirchhoff’s considerations, one may treat vorticity as a sum of point 
vortices and determine the velocity field created by the vortices (this velocity field 
induces the vortex motion) by means of the Biot-Savart law. This way, one may find 
a set of ordinary differential equations (ODEs) for the location of the vortices, that 
describe vortex-vortex interactions [273]. A more subtle consideration in this con¬ 
text is the “screening” effect of the vortex-vortex interaction by the inhomogeneities 
in the density (due to the presence of the external potential). The latter effect has 
been incorporated in some of the above works by an effective renormalization of 
the interaction prefactor, but a more systematic study of this effect is still lacking 
despite the fact that the relevant more complicated dynamical equations can still be 
systematically derived [274]. 

In a 2D (disk-shaped) BEC confined in a parabolic trap, it is then possible 
to employ variational arguments [275] and combine both effects: vortex precession 
(induced by the trap) and vortex-vortex interactions. This way, the effective dynamics 
of a small cluster of interacting vortices (of potentially same or different charges) is 
described by a system of ODEs for the centers of the vortices. The relevant dynamical 
system possesses two integrals of motion (Hamiltonian and angular momentum) and, 
thus, it is completely integrable for vortex dipoles composed by two counter-rotating 
or co-rotating vortices. Importantly, a theoretical description of vortex trajectories 
was found to be in excellent agreement with pertinent experimental findings [77]. 
Notice that apart from vortex dipoles, the same methodology has been used in cases 
of vortex clusters composed by more than two vortices. Clusters involving e.g. 4, 6, 
8 vortices of alternating charges in polygonal form [271] or 4, 5, 6 vortices in a linear 
configuration are currently challenging to produce experimentally (and are found 
to be dynamically unstable when possessing more than 4 vortices in a polygonal 
shape, or 3 or more vortices in a linear configuration [276]). On the other hand, 
producing such clusters with a controllable number of vortices of the same charge 
is straightforward (see the second item in Ref. [77]). Please note that the system 
possesses a number of intriguing symmetry-breaking bifurcations [277], which can 
be explained even in analytical form through a linear stability analysis [278]. 

Apart from single vortices and small vortex clusters, there has been much inter¬ 
est in vortex lattices in rapidly rotating condensates. Such configurations consist of a 
large number of ordered lattices of vortices, arranged in triangular configurations, the 
so-called Abrikosov lattices [279]. The first BEC experiments reported observation 
of vortex lattices consisted of just a few (< 15) vortices [280], but later on it was pos¬ 
sible to nucleate experimentally and maintain vortex lattices with over 100 vortices 
[281, 282]. Subsequent efforts enabled the observation of intriguing phenomena as¬ 
sociated with these vortex lattices, including their collective (so-called Tkachenko) 
oscillations, as well as their structural phase transitions either under multi-component 
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interaction (a transition from hexagonal to square lattice was observed experimen¬ 
tally in [283]) or in the presence of external potentials (a similar transition was the¬ 
oretically reported in the presence of a square optical lattice in [284]). Lastly, it 
is relevant to mention here that such multi-vortex configurations are at the heart of 
ongoing studies of phenomena including vortex turbulence and more generally non¬ 
equilibrium dynamics of atomic condensates [67]. 

Admittedly, the above discussion contains a rather partial perspective of a field 
that has truly boomed in a remarkable number of research threads and has done so via 
a rather unique cross-pollination with other fields of physics that is simply impossible 
to capture within the confines of the present chapter. Nevertheless, we hope to have 
conveyed some of the main areas of the pertinent studies and the ever-expanding 
(in terms of research groups and themes of study) enthusiasm surrounding this area. 
We now turn our attention to the specific contributions associated with this special 
volume, which touch upon many of the above-mentioned topics. 

4. A SYNOPSIS OF THE ARTICLES INCLUDED IN THE PRESENT SPECIAL ISSUE 

As mentioned above, the current theoretical and experimental work on BEC 
and related topics covers a vast research area. Of course, the papers selected for 
this Special Issue cannot survey all aspects of this work. Most of the papers present 
theoretical results, in compliance with the obvious trend that many more original 
theoretical papers on BEC, than experimental ones, appear in the scientific literature. 
Nevertheless, some articles from the Special Issue present experimental results too, 
as briefly recapitulated below. Some papers deal directly with basic aspects of the 
studies of BEC, while others address different but related topics, such as fermion 
quantum gases, few-boson models, etc. The different settings and problems ad¬ 
dressed in the articles may be categorized as more physical or more mathematical 
ones. 

(1.) M. A. Caracanhas, E. A. L. Henn, and V. S. Bagnato, Quantum turbulence 
in trapped BEC: New perspectives for a long lasting problem 

BEC in atomic gases provides the most natural testbed for exploring turbulent 
dynamics of superfluids. This article [285] offers a review of recent experimental 
and theoretical results on this topic, and a discussion on the directions for the further 
development of studies dedicated to quantum-liquid turbulence. 

(2.) A. Vardi, Chaos, ergodization, and thermalization with few-mode Bose- 
Einstein condensates 

This article [286] considers a system with few degrees of freedom, which rep¬ 
resents “small” Bose-Hubbard (BH) models, namely, BH dimers and trimers, the for¬ 
mer one being reduced to a classical kicked top. In the framework of these systems, 



25 


Bose-Einstein condensation: Twenty years after 


29 


the analysis is focused on aspects of classical dynamical chaos in them, including the 
problems of the onset of ergodicity and thermalization. The energy diffusion in the 
systems' phase space is explored by means of the Fokker-Planck equation. 

(3.) R. Radha and R S. Vinayagam, An analytical window into the world of 
ultracold atoms 

The paper [287] addresses BEC models based on GPEs, in terms of the possi¬ 
ble integrability of these models. The approach develops the known method of trans¬ 
forming the standard integrable form of the NLSE into seemingly complex, but still 
integrable ones, by means of explicit transformations of the wave functions and vari¬ 
ables (x,t). In particular, considered are models which may be integrable, while they 
include complex ingredients, such as the time dependence (management) of the scat¬ 
tering length, and a parabolic potential (expulsive, i.e., of the anti-harmonic-oscillator 
type, rather than the trapping one), with a constant or time-dependent strength. In ad¬ 
dition to the single-component models, two-component systems are considered too, 
with both nonlinear and linear couplings between the components. A number of 
exact solutions are found in such models, including bright and dark solitons. 

(4.) A. I. Nicolin, M. C. Raportaru, and A. Balaz, Effective low-dimensional 
polynomial equations for Bose-Einstein condensates 

The article [288] addresses the derivation and analysis of effective equations 
with reduced (ID and 2D) dimensions for prolate and oblate (cigar-shaped and pan- 
cake-shaped) condensates, respectively. The equations with the reduced dimension¬ 
ality are derived from the full 3D GPE, under the condition of strong confinement 
in the transverse direction(s). The effective equations with polynomial nonlinearities 
are derived in this context. 

(5.) V. I. Yukalov and E. P. Yukalova, Statistical models of nonequilibrium 
Bose gases 

The analysis in this article [289] addresses strongly perturbed BEC, i.e., it goes 
far beyond the limits of the near-equilibrium mean-field theory. In particular, this 
paper makes a contact with the considerations presented in the article by M. A. Cara- 
canhas, E. A. L. Henn, and V. S. Bagnato in the same Special Issue [285], as the anal¬ 
ysis develops a description of strongly excited BEC in terms of a statistical model of 
grain turbulence. 

(6.) H.-S. Tao, W. Wu, Y.-H. Chen, and W.-M. Liu, Quantum phase transitions 
of cold atoms in honeycomb optical lattices 

Optical lattices with different geometries (in particular, 2D honeycomb lat¬ 
tices considered in this article) help to support quantum gases (both bosonic and 
fermionic) in highly-correlated states. The article [290] aims to review recent results 
for quantum phase transitions of cold fermionic atoms in these lattices. In that sense, 
it is an essential addition to the collection of topical chapters on the theme of BEC, as 
results are presented for quantum Fermi gases, rather than for bosons. The analysis 
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combines mean-field considerations with quantum Monte-Carlo computations, with 
the aim to calculate various properties of the systems under consideration, such as 
the density of states, the Fermi surface, etc. Also considered in this article are bilayer 
lattices, in addition to the monolayer ones, and effects of the spin-orbit interaction 
between the fermion components. 

(7.) T. He, W. Li, L. Li, J. Liu, and Q. Niu, Stationary solutions for nonlinear 
Schrodinger equation with ring trap and their evolution under the periodic kick force 

This work [291] addresses solutions of the nonlinear Schrodinger equation 
for a periodically kicked quantum rotator. The model may also be relevant for a 
BEC trapped in a toroidal quasi-ID trap, with a periodically applied potential pro¬ 
file. The analysis is focused on the especially interesting cases of quantum anti¬ 
resonance and quantum resonance, using analytical stationary solutions of the non¬ 
linear Schrodinger equation with periodic boundary conditions. 

(8.) D. A. Zezyulin and V. V. Konotop, Stationary vortex flows and macro¬ 
scopic Zeno effect in Bose-Einstein condensates with localized dissipation 

This article [292] addresses an interesting topic of nonlinear BEC in dissipa¬ 
tive media. A specific setting is considered with flow of the superfluid towards the 
central part of the 2D system, where the loss is concentrated. The model does not 
include any explicit gain, but, nevertheless, it gives rise to stationary global patterns, 
including those with embedded vorticity, due to the balance between the influx from 
the reservoir at infinity and the effectively localized dissipation. The solution is in¬ 
terpreted in terms of the Zeno effect in the dissipative BEC. 

(9.) V. Achilleos, D. J. Frantzeskakis, P. G. Kevrekidis, P. Schmelcher, and 

J. Stockhofe, Positive and negative mass solitons in spin-orbit coupled Bose-Einstein 
condensates 

The article [293] addresses the currently hot topic of solitons in the two-com¬ 
ponent BEC realizing the spin-orbit-coupling effect. The analysis is developed for 
the ID geometry, and relies on the reduction of the underlying two-component GPE 
system to a single NLSE, by means of the multiscale-expansion method. In this way, 
apart from the usual positive-mass bright and dark solitons, negative mass structures, 
namely bright (dark) solitons for repulsive (attractive) interactions are predicted as 
well. The analytical predictions are confirmed by numerical simulations. 

(10.) A. I. Yakimenko, S. I. Vilchinskii, Y. M. Bidasyuk, Y. I. Kuriatnikov, 

K. O. Isaieva, and M. Weyrauch, Generation and decay of persistent currents in a 
toroidal Bose-Einstein condensate 

The topic of persistent superfluid flows in toroidal traps is theoretically ad¬ 
dressed in the article, being motivated by recent experimental observations of this 
effect. This article [294] offers a review of theoretical results on this topic, recently 
produced by the present authors. In particular, special attention is paid to the phe¬ 
nomenon of hysteresis in this setting. The analysis is performed by means of numer- 
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ical solutions of the 3D GPE. Some related 2D settings are considered too. 

(11.) M. Galante, G. Mazzarella, and L. Salasnich, Analytical results on quan¬ 
tum correlations of few bosons in a double-well trap 

This work [295] deals not with BEC proper, but rather with sets of few bosons, 
the number of which is N = 2, 3, or 4. The bosons are trapped in a double-well 
potential, which is also used in many experimental and theoretical studies of BEC, 
such as those dealing with bosonic Josephson oscillations. Eventually, the system is 
reduced to the simplest two-site truncation of the Bose-Hubbard model, which has 
something in common with the setting considered in another article included into 
this Special Issue, the one by Vardi [286]. The analysis aims to find exact ground 
states for these few-boson sets, and study variation of their characteristics, such as 
the energy and entanglement entropy, as functions of system's parameters. 

(12.) V. Bolpasi and W. von Klitzing, Adiabatic potentials and atom lasers 
The article [296] addresses the topic of the design of atom-beam lasers. A 
detailed analytic model of the trap is presented and the flux of the atom laser is 
determined. The analytical results are found to be in good agreement with recent 
experimental data. The analysis is focused on the harmonic-oscillator trapping po¬ 
tential for the BEC, from which the laser beams are emitted. Gravity is taken into 
regard too. 
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